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The vast majority of mutations are deleterious, and are eliminated 
by purifying selection. Yet in finite asexual populations, purifying 
selection cannot completely prevent the accumulation of delete- 
rious mutations due to Muller's ratchet: once lost by stochastic 
drift, the most-fit class of genotypes is lost forever. If deleterious 
mutations are weakly selected, Muller's ratchet turns into a mu- 
tational "meltdown" leading to a rapid degradation of population 
fitness. Evidently, the long term stability of an asexual popula- 
tion requires an influx of beneficial mutations that continuously 
compensate for the accumulation of the weakly deleterious ones. 
Here we propose that the stable evolutionary state of a population 
in a static environment is a dynamic mutation-selection balance, 
where accumulation of deleterious mutations is on average off- 
set by the influx of beneficial mutations. We argue that this state 
exists for any population size N and mutation rate U. Assuming 
that beneficial and deleterious mutations have the same fitness 
effect s, we calculate the fraction of beneficial mutations, €, that 
maintains the balanced state. We find that a surprisingly low e 
suffices to maintain stability, even in small populations in the face 
of high mutation rates and weak selection. This may explain the 
maintenance of mitochondria and other asexual genomes, and 
has implications for the expected statistics of genetic diversity in 
these populations. 

asexual evolution | Muller's ratchet | adaptation 

Purifying selection maintains well-adapted genotypes in the 
face of deleterious mutations [1]. Yet in asexual popula- 
tions, random genetic drift in the most-fit class of individu- 
als will occasionally lead to its irreversible extinction, a pro- 
cess known as Muller's ratchet [2, 3]. The repetitive action 
of the ratchet leads to the accumulation of deleterious muta- 
tions, despite the action of purifying selection. This ratchet 
effect has been extensively analyzed [4, 5, 6, 7, 8] and has 
been observed in experiments [9, 10, 11, 12] and in nature 
[13, 14, 15]. In small populations when deleterious mutation 
rates are high or selection pressures are weak, the ratchet can 
proceed quickly, causing rapid degradation of asexual genomes 
and a collapse in population size: the so-called "mutational 
meltdown" [16, 17, 18]. Hence, Muller's ratchet has been de- 
scribed as a central problem for the maintenance of asexual 
populations such as mitochondria [19]. Avoiding this muta- 
tional catastrophe is thought to be a major benefit of sex and 
recombination (see [20, 21] for reviews). 

In the absence of recombination [22] or epistasis [23], the 
only forces that can check the accumulation of deleterious mu- 
tations are back and compensatory mutations. These mu- 
tations are typically assumed to be rare and are either ne- 
glected in models of the ratchet, or assumed to only slightly 
slow its rate [1, 24, 25]. However, extensive experimental 
work [26, 27, 28, 29, 30, 31, 32] and theoretical expectations 
[33, 34, 35, 36] suggest that the number of available back and 
compensatory mutations increases as a population accumu- 
lates deleterious mutations and declines in fitness. The in- 
crease in beneficial mutation rates with declining fitness will 
eventually stop the progress of Muller's ratchet, as illustrated 
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Fig. 1. The relationship between the fraction of mutations that are beneficial, 
e, and the absolute fitness of the population, a). A poorly adapted population 
(green) with o) < a>c and e > £c, will adapt towards the dynamic equilibrium 
"Attractor" state (gray) with higher fitness and lower e. Conversely an "over- 
adapted" population (red) with a> > a>c and e < £c, will decline in fitness towards 
the dynamic equilibrium state due to Muller's ratchet. 



in Fig. 1. Furthermore, since back and compensatory muta- 
tions are selectively favored, they can balance deleterious mu- 
tations even while they are relatively rare, thereby maintaining 
a well-adapted population. 

Despite the potential relevance of compensatory mutations, 
only a few studies have quantitatively analyzed the simultane- 
ous accumulation of both beneficial and deleterious mutations. 
These include simulations by Antezana and Hudson [37] and 
Wagner and Gabriel [38], experimental work by Silander et al. 
[28], and both simulations and analytical work by Rouzine et 
al. [39, 40]. In each of these analyses, simulations or exper- 
iments showed that the mean population fitness settles at a 
steady-state value co^ — the point where beneficial mutations 
arise often enough to compensate for newly-fixed deleterious 
mutations. 

While the absolute fitness, cDc, of the population in this dy- 
namic equilibrium state is model dependent, the existence of 
the stable steady state is generic. Intuitively, suppose a fraction 
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e of mutations are beneficial, and the critical fraction needed 
to maintain stability is e^. In a poorly adapted population, a 
higher fraction of mutations will be beneficial, e > e^. This ex- 
cess of beneficial mutations will push the population towards 
higher fitness. At the same time, adaptation will deplete the 
available pool of beneficial mutations, until is reached, as 
shown in Fig. 1. Conversely, in an "over-adapted" population, 
we expect e < e^, so that deleterious mutations dominate, re- 
ducing population fitness until is recovered. Thus we expect 
the dynamic mutation-selection balance point to be a stable 
evolutionary "attractor." In contrast to oJc, the critical pro- 
vides a general measure of how well adapted a population is, 
independent of the model of compensatory mutations. 

Below, we show analytically how depends on population 
size, mutation rate, and strength of selection and compare our 
predictions to numerical simulations. We show that even for 
high mutation rate and small population size, the stable state 
occurs at high fitness, where most mutations are deleterious. 
MuUer's ratchet notwithstanding, selection enables rare ben- 
eficial mutations to compensate for more frequent deleterious 
mutations, and to maintain a well adapted population. 



Model 

Our analysis is based on the standard discrete generation 
Wright-Fisher model with population size fixed at N. We 
assume that beneficial and deleterious mutations increase or 
decrease fitness by s (where s > by convention). Each indi- 
vidual has some number, fc, deleterious mutations relative to 
the perfectly adapted state. We define fc to be the mean number 
of deleterious mutations per individual, so that an individual 
with /c deleterious mutations has fitness s(fc - fc) relative to the 
mean, as shown in Fig. 2. The population can be character- 
ized by the number of individuals in each of these fitness 
classes. We define \1 to be the total mutation rate, and e to be 
the fraction of mutations that are beneficial, so that the bene- 
ficial mutation rate is (Jf, = (Je, and the deleterious mutation 
rate is (J,( = - e). 

Although mutations and genetic drift are stochastic pro- 
cesses, in the bulk of the distribution — where fitness classes 
contain many individuals — the population dynamics are well- 
captured by a deterministic approximation. On average, the 
number of individuals in fitness class fc evolves as: 



A 



= s(k - k)nk - Urik + Udrik-i + Utrik+i . 



[1] 



The terms on the right hand side of this equation describe the 
different processes acting on fitness class k: (I) the effect of se- 
lection relative to the mean fitness, (II) mutational load, (III) 
deleterious mutations from more-fit individuals, and (IV) ben- 
eficial mutations from less-fit individuals. We assume that the 
selective effect of a single mutation is small, s ^ 1. 

In general, we expect e to increase with k, since less-fit 
genotypes have a higher frequency of beneficial back and com- 
pensatory mutations. However, at any given time there will 
only be a relatively narrow range of k present in the popula- 
tion. Thus in Eq. (1) we neglect the fc-dependence of e; we have 
tested this approximation against a simulation incorporating 
the Jc-dependence of e (see SI). We use this approximation to 
compute the £c and corresponding distribution n/c at which the 
population stays on average at a constant fitness. As we have 
argued above, the stability of this equilibrium state is ensured 
by the fc-dependence of e on larger scales. 
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Fig. 2. (A) Schematic illustration of the fitness distribution within a population. 
We refer to the most-fit class as the "nose" of the distribution. (B) A typical re- 
alization of the stochastic dynamics at the nose. In the event of extinction of the 
nose, «! becomes the new nose and relative nose fitness decreases by s. Con- 
versely, when a more-fit nose is established the relative nose fitness increases 
by s. 



Results 

In steady state, the fitness distribution stays constant on av- 
erage: ^ = Oforallfc. Solving Eq. (1) for the steady state in the 
case of e = leads to the familiar mutation-selection balance 
[1], with a Poisson distribution of fitness: = Ne"^ ^, where 
A = Lf/s and A: = corresponds to the most-fit class. 

This deterministic approximation corresponds to an infinite 
population where MuUer 's ratchet (i.e. the stochastic extinction 
of the fittest genotype) does not operate. In the absence of such 
stochastic extinction, any beneficial mutation from the fittest 
class would move the distribution towards higher fitness, in 
conflict with the steady state assumption. Hence, in infinite 
populations, steady state can only be achieved with = 0. 
Conversely, in a finite population where the fittest class can be 
lost due to genetic drift, MuUer's ratchet will eventually lead 
to a decrease in fitness if no beneficial mutations are available. 
Thus we must have > to be in the dynamic steady state. 

A correct description of the dynamic equilibrium state 
therefore requires a suitable treatment of genetic drift. These 
stochastic effects are particularly important in the most-fit edge 
of the distribution, which we call the "nose," where the number 
of individuals with a particular fitness is small. Our analysis is 
based on matching a stochastic treatment of the nose with a de- 
terministic description of the bulk of the fitness distribution (as 
governed by Eq. (1) with e > 0). The stochastic dynamics of the 
nose are determined by two competing processes: the random 
extinction of the most-fit class, versus the establishment of a 
new more-fit class due to a beneficial mutation. At stationarity, 
the rates of these two processes have to be equal; this condi- 
tion determines the number of individuals in the nose class. 
We match this stochastic condition with the number of indi- 
viduals in the nose class calculated based on the deterministic 
distribution. This determines the critical fraction of beneficial 
mutations £c as a function of the population parameters. 
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Fig. 3. Comparison of the fitness distribution obtained in simulations with the 
analytic result Eq. (2) and a Gaussian with equal mean and variance. The sim- 
ulation result is the median fitness distribution observed over 10* generations. 



The shape of the fitness distribution. We begin with a determin- 
istic analysis of the shape of the fitness distribution at steady 
state. Since only relative fitnesses matter, it is convenient to set 
the origin of k such that the mean of the fitness distribution is 
k = A. We solve for the steady state of Eq. (1) using Fourier 
analysis (see SI). We find 

= Ne*'^(i-2f)-*:iog ^/WJ|^(2a), [2] 

where }k denotes the Bessel function of order k and we have 
defined a = A -y/e(l - e). A comparison of this solution to sim- 
ulation results is shown in Fig. 3. For k in the vicinity of 
k, this distribution is approximately Gaussian with variance 

= A(l - 2e). While the the low fitness tail of the distribu- 
tion decays less rapidly than a Gaussian, the high fitness side 
decays more rapidly than a Gaussian. 

It is important to note that in the deterministic limit there 
is no true stationary state for arbitrary e > 0; this manifests 
itself in the loss of positivity of the solution given by Eq. (2) for 
high-fitness classes above the position of the nose class. While 
the position of the nose, which we will call k^,, requires a more 
careful treatment, Eq. (2) gives a very accurate description of 
the bulk of the fitness distribution (i.e. for k > ki,). 

To determine the position of the nose, fc*, we observe 
that if classes with < fc» carry no individuals, fitness class 
fc* does not receive an influx of mutations from the more-fit 
class ki, - 1. That is, term III in Eq. (1) is absent, yielding 
{k-ki, - A)nki, + AeH/c,+i = 0. This along with Eq. (2), gives the 
following equation for fc* 

fc*/,,(2a)= VS/t,+i(2a) [3] 

This must in general be determined numerically. However, in 
SI we show that 

(A^e foreA^ « 1 

* ~ \2A Ve(l - e) for eA^ » 1. 

To gain some intuition into the significance of fc*, note that 
A:* = in an infinite population with = 0. In this case, the 
fitness distribution is the familiar Poisson distribution and the 
fittest class contains a fraction e"'^ of all individuals. In a finite 
population, Ne^'^ will be smaller than one for sufficiently large 
A, which implies that at steady state the mutation free geno- 
type is typically absent and therefore fc* is larger than 0. This 
adjustment of fc* for finite population size is reflected in Eq. (4), 
which states that the most-fit class gets closer to the population 
mean as e increases. 



The stochastic matching condition. We seek to determine the 
£ = £c at which the distribution is stationary for a given finite 
N. Since the most-fit class is populated by a comparatively 
small number of individuals, fluctuations due to genetic drift 
may change its occupancy significantly. In particular, the pop- 
ulation at the nose can go extinct, in which case the next fitness 
class becomes the new nose. Alternatively, a lucky beneficial 
mutation can cause the nose to advance by one class. Station- 
arity can therefore only be achieved in an average sense: the 
rate of advancing the nose has to equal the rate of extinction of 
the nose. 

Since these two processes depend sensitively on the num- 
ber of individuals in the nose, requiring the equality of extinc- 
tion and establishment rates determines the average size of the 
population of the nose, Wt^, as a function of U, s, and e. We 
match this n^^ to the deterministic solution determined above 
tofindec(N,s, U). 

The nature of the dynamics at the nose depends qualita- 
tively on W)c^ . If is large enough that its dynamics are domi- 
nated by selection, it is rarely lost and MuUer's ratchet is slow. 
Conversely, if w^^ is small, the class turns over neutrally and 
is easily lost and reseeded. We begin by considering the slow 
ratchet regime where M/c, is relatively large, and then turn to 
the opposite fast ratchet regime. 

Slow ratchet regime. In the regime where rit^s > 1, the fittest 
class is only rarely lost due to drift. Note that this regime cor- 
responds to eA^ < 1, which implies nt^ ~ Ne"'^ and fc» ~ 0. 
In this regime the mean extinction rate for the nose, r_, is 
due to rare large fluctuations which overcome the "restor- 
ing force" due to selection trying to preserve mutation selec- 
tion balance. Estimates of this extinction rate have been ob- 
tained via diffusion theory [1, 6, 4] and have the form (see SI) 
r_ a; e")'^"*. ys -yjysni^^ In, where y is a parameter characterizing 
the effective strength of selection on the fittest class, introduced 
by Haigh[l]. 

To impose our stochastic condition, we must also calculate 
the rate at which a new more-fit class establishes, r+. This is 
given by the product of the rate at which new beneficial mutants 
are generated and the probability they establish, r+ = rii^^ UtPcst- 
Since fc« » in this regime, we have Pest ~ 2s. Equating the two 
rates yields the stochastic condition 

ec = |j . , :. [5] 

^ysNe-^n 

This condition suggests that in the slow ratchet regime 
is a function of a combination of population parameters: 
A£c{U,'N,s)ly^ ~ f{ysNe~'''). We find that contrary to earlier 
studies [1, 4] constant y is not consistent with simulations. In- 
stead, as shown in Fig. 4A the simulation data is well explained 

hyy = 1/ Va. Our empirical A-dependence is however consis- 
tent with the approximation y ~ 0.6 used by Gordo et al. [4], 

within the range of parameters (0.4 < 1/ Va < 0.8) addressed 
by their study. Since we consider a much broader range of 

parameters (0.1 < 1/ Va < 1), accounting for the dependence 
of 7 on A = Lf/s becomes important. Our analytical solution 
Eq. (5), shown as dotted line in Fig. 4A, is in good agreement 
with the simulation data across this wide range of parameters. 

In this regime, depends exponentially on Nse~'\ and 
hence rapidly approaches zero as A becomes small, approach- 
ing the infinite population limit where MuUer's ratchet does 
not operate. 

Fast ratchet regime. For larger mutation rate, smaller pop- 
ulation size, or weaker selection, the occupancy of the most 
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Fig. 4. The critical proportion of beneficial mutations, £c, across both slow and 
fast ratchet regimes for many combinations of U, s and N for the two regimes. 
Simulation results are shown as points, with color Indicating A = ^. (A) Slow 

ratchet regime. We show A^e^ as a function of -^sNe~'^ for A^e^ < 1- The 

dashed black line Is the analytic solution given by Eq. (5). (B) Fast ratchet 



regime. Note that £c is a function of 



logNs 
A 



and even In small populations only 



approaches ec ~ j for very large A. The dashed black line shows the asymp- 
totic solution £c in the large population size limit (ejA^ » 1), Eq. (7). The inset 
compares the numerical solution from Eq. (6) to simulation data. 



This condition, together with the solution Eq. (2) for the distri- 
bution, allows us to determine £c(,N,s, U). As shown in Fig. 4, 
our solution Eq. (6) is in excellent agreement with simulations 
for most of the data. However, our solution overestimates 
for populations with small population size and small A that still 
satisfy n^^s < 1. Abetter description in this regime will require 
a more careful analysis of fluctuations beyond the nose. 

Large population size limit. Large populations can maintain 
a well adapted genome with « | even for moderately large 
A. In the limit of large populations and large A, with the limit 
taken so that £( A^ » 1 while Cc \, the matching condition 

reduces to yi^logec = A'Mog^Nse"'^*^"^*^''). This equation has 
the approximate solution 

log Ns 



log^[2/log(z-i)] 



22 = 1- 



A 



[7] 



This result can be made more precise through iteration (see SI 
for derivations and comparisons to simulation data). Note here 
that Ec depends only very weakly on Ns and on A, in contrast to 
the slow ratchet regime where £c declines exponentially with 

High mutation rate or weak purifying selection. In this limit, 
the dynamic equilibrium tends towards the state where ben- 
eficial and deleterious mutations are equally frequent, £ = \- 
This corresponds to a population that can no longer maintain 
a well-adapted state. In the SI, we derive an approximate ex- 
pression for in the fast-ratchet regime in the limit of large 
A: 



/ 3 \^'^ 

(iAi^sH • 



[8] 



^'^2 liA 

From this expression we can immediately determine the point 
at which the population can no longer maintain a well-adapted 
state, as defined by a maximum fraction of beneficial mutations, 
Gmax- In order to have £c < e„ax, we require 



A 



logNs 



0.5- 



[9] 



valid for e^ax ^ 0.2. Note that for e^ax = 1/4, the RHS of this 
expression is 48; for ejnax = 1/3, the RHS is 162. Thus even small 
populations with Ns > 1 can maintain relatively well-adapted 
genomes in the face of high mutation rates and weak selection 
(i.e. when U » s). 



fit class decreases, thereby increasing its rate of extinction. 
Consequently, a higher rate of beneficial mutations (larger e) 
is required to match the extinction rate. The resulting rapid 
turnover of the population at the nose leads to the failure of the 
quasi-static approximation we used in the slow ratchet regime. 

As the occupancy of the nose decreases, in particular when 
nkt.s < 1, the dynamics of the fittest class is governed by drift. 
The rate of extinction, r., can therefore be estimated from neu- 
tral diffusion: r_ ~ l/{2nk_^). The rate at which a new more-fit 
class is established is given by the same formula as before, 
r+ = £Unk_^Pfsf However, now refers to the probability 
that a new mutant lineage reaches 2(ki!+i)s individuals. At this 
point, the lineage crosses over from stochastic to deterministic 
dynamics, entering the domain described by the deterministic 
solution Eq. (2). Thus P^f = 2(fc» -I- l)s. Equating = r+ yields 
the stochastic condition: 



Discussion 

In infinite populations it has long been recognized that the 
balance between mutational pressure and purifying selection 
leads to a fitness equilibrium [1, 41]. Our analysis demonstrates 
that such an equilibrium also exists in a finite population, de- 
spite the action of genetic drift and MuUer's ratchet. Previous 
work on these finite populations has focused primarily on the 
rate of MuUer's ratchet in the absence of beneficial mutations, 
and has generally assumed that back and compensatory muta- 
tions are rare enough to be neglected [42, 43, 6, 4, 5, 25]. How- 
ever, we expect that beneficial mutations will become more 
common as a population accumulates deleterious mutations. 
Our analysis has illustrated how these beneficial mutations 
balance MuUer's ratchet in asexual populations, to create a dy- 
namic equilibrium state. 

We have shown that for any population size, mutation rate, 
and selection pressure, this dynamic equilibrium state is char- 
acterized by the proportion of beneficial mutations, £c, neces- 
sary to balance the ratchet. This is a measure of how optimal 
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Fig. 5. The fraction of beneficial mutations necessary to maintain the dynamic 
mutation selection balance in a population with parameters Ns and A = j. The 
dotted line separates the slow and fast ratchet regimes. Experiments with various 
model organisms with different population parameters are represented as points. 



a genome a finite population can maintain in the presence of 
mutations, selection, and genetic drift. Larger e^. indicates that 
the population needs more beneficial mutations to be in the 
stable state, and hence must be less fit, and vice versa. In the 
limit £c —> 0, the population is perfectly adapted, and no bene- 
ficial mutations are possible. The opposite case of mutational 
meltdown, where a population cannot maintain an adapted 
genome, corresponds to — > 1/2, where a random mutation 
is as likely to be beneficial as deleterious. There are two qual- 
itatively different regimes for e^/ as shown in Fig. 5. In the 
slow ratchet regime, selection stabilizes the nose, so is small 
and depends exponentially on population parameters Ns and 
A. On the other hand, in the fast ratchet regime, must be 
larger to balance the faster ratchet, and depends more weakly 
on the population parameters Ns and A. The boundary of the 
two regimes is roughly given by Nse~^ = 1. 

We note that Rouzine et al. [39, 40] studied population dy- 
namics under the combined action of purifying and positive 
selection in a model similar to ours, finding both simulation 
and analytical evidence for the existence of the dynamic bal- 
ance state. However, their analysis focuses on the rate of the 
ratchet (with e <K e^) or the rate of adaptation (with e » e^), 
and does not apply as e — > e^. Our work, by contrast, focuses 
exclusively on the dynamic equilibrium state and the critical £c 
required to maintain it. 

Our analytic results for in the dynamic equilibrium were 
made possible by the explicit calculation of the shape of this 
fitness distribution, given in Eq. (2) and illustrated in Fig. 3. 
With increasing U/ s and decreasing N this distribution deviates 
from the Poisson distribution for the classic mutation-selection 
balance [1], most notably reducing the fitness of the top class, 
k^,, relative to the mean. Examining the properties of the dy- 
namic mutation-selection equilibrium over the full parameter 
range, shown in Fig. 5, has revealed a strong asymmetry be- 
tween beneficial and deleterious mutations. For example, for 
a small population Ns ~ 10^ with high mutation rate U/s ~ 8, 
just 2% beneficial mutations are enough to counteract the effect 
of deleterious mutations. This indicates that purifying selec- 
tion is remarkably effective even for conditions where MuUer's 
ratchet would proceed extremely quickly in the absence of back 
and compensatory mutations. 

Our analysis implies that populations with e > should 
adapt, while populations with e < should decline in fitness. 
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Experimental evolution of model organisms in controlled lab- 
oratory environments appears to be consistent with this expec- 
tation. In particular, Silander et al. [28] showed that bacterio- 
phage (pX174: converged to a population-size dependent fitness 
plateau, as our model would predict. For their population pa- 
rameters (s ~ 0.08, U ~ 0.13) our theory predicts ranging 
from 2% to 0.2% for their experiments with N ranging from 
100 to 200, consistent with the beneficial mutation rates they 
infer. Their experiments with lower population sizes (N < 30) 
are in the Fast ratchet regime and estimating requires further 
analysis of fluctuations beyond the nose. 

Other experiments with vesicular stomatitis virus found 
that large populations adapt while small populations melt, also 
consistent with our analysis [44]. Similar results have been ob- 
served in yeast [45, 46]. In the long-term evolution experiments 
of Lenski and collaborators in E. coli, parameters are such that 
we expect the critical Cc to be very small, ^ 10"^*^, consis- 
tent with their observation of continuous adaptation even after 
tens of thousands of generations. On the other hand, mutator 
strains of E. coli shows decrease in fitness [47], which could 
be used to further test the model. In Fig. 5, we show rough 
estimates for the parameter regimes in which each of these 
experimental systems lie. 

Population parameters for HIV place it close to the bound- 
ary between the slow and fast ratchet regimes (Fig. 5), so that 
its evolutionary dynamics depend sensitively on the strength of 
selection. This intriguing case highlights two important issues. 
First, our analysis has assumed that all mutations have the same 
selective effect s, illustrating how well a particular strength 
of purifying selection can maintain well-adapted states. In 
general, however, mutations have a range of selective effects. 
Mutations with a particular characteristic strength may tend 
to dominate the dynamics, if more strongly selected mutations 
are quickly eliminated by selection while more weakly selected 
mutations are effectively neutral given interference from selec- 
tion on the other mutations. However, the effect is non-trivial 
and generalizing our results to an arbitrary distribution of mu- 
tational effects remains an important avenue for future work. 
Second, our analysis has focused on asexual populations, and 
neglected recombination. Since even weak recombination has 
the potential to significantly slow MuUer's ratchet, it would be 
interesting to generalize our dynamic balance to include its ef- 
fects. Increasing recombination rates would presumably allow 
the population to maintain better-adapted genotypes in this 
steady state. 

We have argued that the dynamic mutation selection bal- 
ance is a stable evolutionary attractor. If the population size 
is independent of e, this is immediately obvious from Fig. 1. 
Under some conditions, however, the population size might be 
a function of absolute fitness and hence of e. MuUer's ratchet 
might then lead to a decrease in fitness, resulting in a smaller 
population size, which in turn leads to a faster ratchet. This 
phenomenon of mutational meltdown might cause rapid pop- 
ulation extinction [16]. Eq. (8), however, shows that the degree 
of adaptation depends only weakly on the population size and 
a stable equilibrium can be maintained even if N depends on e 
rather strongly. 

The dynamic balance state has important implications for 
patterns of molecular evolution. Recent analysis of the effects 
of purifying selection on the structure of genealogies has sug- 
gested that MuUer's ratchet plays a crucial role in determin- 
ing the structure of genetic variation in asexual populations 
or on short distance scales in the genomes of sexual organ- 
isms [48, 49]. However, these expectations must be revised if 
populations exist instead in the dynamic state in which both 
beneficial and deleterious mutations fix, without any contin- 
uous net accumulation of deleterious mutations. This means 
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that even though no change in fitness occurs, signatures of both 
positive and negative selection are likely to be found in patterns 
of molecular evolution, as has been suggested by earlier studies 
[50, 37]. We argue that this state is the natural null expectation 
for the effects of mutations and purifying selection on patterns 
of genetic variation; efforts to look for positive selection which 
represent "true" adaptation should look for deviations from 
this situation. 

Finally, we note that although the dynamic balance we have 
analyzed is stable, a population will typically fluctuate aroimd 
this steady state. A few beneficial mutations may become es- 
tablished at the nose by chance, leading to a temporary increase 
in mean fitness which is later balanced by a reduction in e at 
this higher fitness, restoring the population to the equilibrium 
state. Conversely, a few clicks of MuUer's ratchet will occasion- 
ally lead to temporary reductions in population fitness before 



the corresponding increase in e restores the steady state. These 
fluctuations aroimd the dynamic balance may be important 
to the clonal structure of the population, and hence are likely 

to play a key role in patterns of molecular evolution and in 
understanding the effects of recombination. 
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Supplemental Information 

Simulation methods. Our simulations are done using a cus- 
tom written Python code available on request. We implement 
a discrete time Wright-Fisher model where the population is 
represented by a vector with elements corresponding to the 
number of individuals in fitness class k. Each generation con- 
sists of separate selection and mutation steps. To implement 
selection, the vector n^t is multiplied by e***^"^'"" to obtain a vec- 
tor nt containing the expected number of offspring in class k. 
Here (k - k) is the fitness of class k relative to the population 

mean, and a = ^"^"'^ maintains an approximately constant 
population size aroimd N. 

Our model is simplified by using beneficial and deleterious 
mutation with effect ±s and mutation rate lie and Lf(l - e), re- 
spectively. To implement mutation, we calculate the probabil- 
ity P{i, j) of a genome being hit by ; beneficial and /' deleterious 
mutations, which are Poisson distributed. This mutation ma- 
trix is then applied to hk, i.e. the parts of fik are moved up or 
down according to the net number of mutations they accrued 
in this time step. Having constructed the expected number 
of individuals in fitness class k after selection and mutation, 
we draw a population sample from each class from a Poisson 
distribution. If necessary, tiie population is re-centered in the 
discrete vector tik- This prevents occupied classes from rim- 
ning off the grid due to accimiulated increase or decrease in 
their absolute fitness. 

The above process is repeated for a specified number of 
generations. The speed of adaptation and other features of the 
dynamics are measured after an equilibration time to remove 
transient effects from the initial conditions. In the parameter 
regimes studied, we found that 10^ generations was generally 
sufficient. 

To solve for edU, N, s), we re-run the simulation while iter- 
atively adjusting e to get v as close to zero as possible. Since v 
is a stochastic quantity, Cc can only be determined with limited 
accuracy, but this can be improved by increasing the number 
of generations for each run. We have also run simulations in 
which e varies with fitness (i.e. with k) to demonstrate that the 
population indeed evolves towards the values of k at which 
they have the appropriate Cc- 

Stationary distribution in the presence of deleterious and 
beneficial mutations. The steady state distribution satisfies 

= (k- k)rn: - Xrik + A(l - e)nk+i + AeTit-i [10] 

which we solve by introducing a Fourier transform 



k 



which yields an ODE in h{q) 

d 



dq 



n(q) = - e)e-"' + ee"f]n(q). 



[11] 



[12] 



The ODE is readily solved to give 

In n(<?) = A [(1 - - ee"'] - A(l - 2e) [ 13 ] 

where the last term is a constant added to enforce normaliza- 
tion E^nt = 1- 

Fourier transforming back to solve for w^^ gives 



This can be re-written by shifting the integration contour 

upward into the complex plane by Q, = In - e)/e, i.e. 
change integration variable q = q' - iQ, which yields 



m = e 



-A(l-2e)-*:log ■ 



[15] 



^ g-A(l-2e)-HogV¥;^(2a) 



where a = \ ^|e{\ - e) and the Bessel fimction is identified from 
its integral representation [51]. 

While this solution holds without any restriction on param- 
eters, it is instructive to examine the e — > Umit, for which we 
have a regular perturbative expansion in eA^ « 1: 



„ _ p-A(l-2e)-J:log[A(l-e)] y » (-l)'"^' 
A2e(l-e) ^ 



~ k\ ' ^ 



[16] 
[17] 



which in the e — > recovers the familiar Poisson muta- 
tion/selection balance. 

Rate of extinction of the nose in slow ratchet regime. The 

rate of extinction for the nose can be calculated using the 
stochastic dynamics of the fittest class, which can approxi- 
mately be modeled by a diffusion equation: 

difix, t) = -d, [Di{x)p{x, t)] + d\ [D2(x)p(x, i)\ [18] 

where Tp(x) = p(nk^,/N = x) is the probability distribution of 
the nose occupancy and Xk^, = e''^^" is the mean fraction of the 
population, Di{x) = sx{l - x/ki,), and Dzix) = xiX - x)/2N. Both 
X = 1 (fixation) and x = (extinction) are absorbing boimdary 
conditions [6]. The density function for fixation/extinction time 
(p(t; x) satisfies the backward Koknogorov equation 



dt(p{t;x) = Di{x)dx(t>{t;x) + D2ix)dl(p{t;x) 



[19] 



where t is the time interval between initial state x and fixation 
or extiaction. The mean time to fixation/extinction starting at 
X = y at t = is given by t{y) = t(p{t; y)dt, which satisfies 

-1 = D,{y)dyt(y) + D^iy^Ky). [20] 

, yields for the 



Using the integrating factor (p{x) = e 
mean time to extinction of the fittest class: 



. Jo D2{z) 



di 



[21] 



The second integral in the limit Ns » 1 can be approximated 
as 



2N r dC 



(1-Q e'"* 
OW) 



■ 2N f dC- 



[22] 

Using a = Nsxk^,,nk^, = Nx^^, and changing the variables, 
<#> = y/xk^, - 1,2 = C,/Xk^, - 1 yields: 



t{Xk,tar) = 2n, 



dcpe"^" 
1 J(j) 



l/kt-l 



dz- 



[23] 



For a » 1 the integral the second integral can be readily ap- 
proximated yielding 



nk = e 



-A(l-2e) 



J ^ ^''P + ^(1 - - ^ee"'] • [ 14 ] ffe* ) ~ 2nfc, J d(pe"'l^ [erf( ^^a^) - erf( ^/a(p)\ [ 24 ] 
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where p = — 1- Using (p^ = 1- 6^ la allows the evaluation of 
the integral in the a » 1 limit and yields 



Assuming an exponential distribution for time of extinction 
of the fittest class, the rate of extinction is given by r_ = l/t{xk_^ ), 
which is used in the Main Text. 



Asymptotic expressions for The matching condition used 
to determine in the fast-ratchet regime can be rewritten as 
follows. 

1 :Ne-'^<i-2'^'+'^'°s^/,,(2a). [26] 



In the limit of a » 1, the zero of the Bessel function is approxi- 
mately at fc* + 1 ~ a and the Bessel fimction at fc* evaluates to 
roughly ~ k^^ (Abramowitz 9.3.33 and 9.1.27). Hence we have 



A(l - 2e) - A Ve(l-e)log - 



logNs 



[27] 



where we have neglected powers to the ^ and (3(1) factors 
inside the logarithm. In the limit e « 1 (but A^e » 1), the 
matching simplifies to 

1 -2e+ Veloge = A'MogNs [28] 

which simplifies further to 

Veloge ~ A"^log(Ns)- 1 = -2z [29] 

where we defined z for convenience. This can be solved for e: 

z^ 

[30] 



W(-z)2 (log(z)-log(-log(z)))2 

where W(x) is the -1 branch of Lambert's W-function, i.e. the 
solutions of W(x)e'^^^^ = x. The linear correction in Eq. 28 can 
be incorporated iteratively. 



e/+i 



(2 + ejf 
W(-z + e,)2 



This iteration converges for small C and e. At larger C, the 
branch of the W{x) function is lost. The result £2 obtained after 



the first iteration, starting with Eq. 30, is also shown in Fig. 6 
as dashed red line. 

The other limit that is amenable to analytic calculations is 



[25] the limit e —> 1/2. To this end, we define 6 - 



the right hand side of Eq. 28 

l-2e- Ve(l-e)log^'^ ^ 
From this, we find 



- £ and expand 



44, 



36' + ^6'+0(6^) = A-i logNs. 

[32] 



^-(^log(Ns) 



1/3 



[33] 



This expression is compared to simulation results in Fig. 6. The 
term proportional to 5^ can again be included by iteration and 
the result is shown as the dashed green line in Fig. 6. 
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Fig. 6. Dynamic mutation-selection balance. The figure shows £c for many com- 
binations of U, s and JV as a function of fj log Ns, while the color codes for 
[31] ectfs"^. If £cU^s~^ a> 1, is solely a function of jj logNs and is well de- 
scribed by the numerical solution of Eq. 6 shown as a black line. The asymptotic 
approximations for large £c (Eq. 33) and small £c (Eq. 30) are shown as green 
and red lines, respectively The dashed lines correspond to the more accurate 
version mentioned in the main text. 
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